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TECHNICAL NOTE 3722 

GENERAL THEORY OF WAVE-DRAG REDUCTION FOR COMBINATIONS 
EMPLOYING QUASI-CYILNDRICAL BODIES WITH AN APPLI- 
CATION TO SWEPT-WING AND BODY COMBINATIONS 1 
By Jack N. Nielsen and William C. Pitts 

SUMMARY 


The wing-body interference theory of NACA TN 2677 applied to symmet- 
rical' 'wings in combination with quasi-cylindrical bodies permits the 
direct calculation of pres sure -distribution changes produced by body 
shape changes. This theory is used to determine the relative magnitude 
of the wave-drag reduction produced by changes in' cylinder cross-sectional 
area and that produced by changes in cross-sectional shape (without change 
in area) . The body distortion is expressed as a Fourier series, and an 
integral equation is derived for the body shape for minimum drag for each 
Fourier component. Thus the wave-drag reductions for the various Fourier 
harmonics are independent and additive. 

To demonstrate the use of the method, it is applied to a body with 
a 4-5° swept wing having a chord of three body radii, a span of six body 
radii, sonic leading and trailing edges, and a biconvex airfoil section. 
Calculations of the various drag components for this configuration show 
that the axially symmetric harmonic of body distortion caused a reduction 
of 22 percent of the wave drag of the wing panels on the cylindrical 
undistorted body. The second harmonic, which causes no change in volume, 
produced a drag saving of 42 percent. The wave-drag saving was obtained 
by decreasing the wing-panel drag by twice as much as the body drag was 
increased. The wing panels have negative wave drag for two harmonics. 
These reductions are accompanied by pressure gradients that are adverse 
for the boundary layer. With regard to body shape, volume is removed 
from the sides of the fuselage and added to the top and bottom. A net 
volume less than the wing volume is removed if the distortions are 
limited to the interval of the root chord. Although the second harmonic 
removes no volume, it causes indentations about five times as deep as 
those due to the axially symmetric harmonic. As a result the wing thick- 
ness for which the full theoretical distortions can be utilized is about 5 
percent or less. 

A very simple theory for the axially symmetric harmonic can be 
derived by neglecting three-dimensional effects. This theory is valid 

^Supersedes recently declassified NACA RM A55B07 by Jack N. Nielsen, 


1955 
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for small chord lengths . An approximate treatment of the three-dimensional 
effects yields a fairly simple theory applicable to longer chord lengths. 
The effects of semispan-radius ratio and chord-radius ratio on indentation 
shape, indentation volume, ‘and drag saving have been determined using- the 
approximate theories for the axially symmetric harmonic. 

An alternate method is developed for determining the body shapes for 
minimum drag for any Fourier harmonic. By this method it is possible to 
obtain for the particular case considered an analytical expression for the 
shape of the axially symmetric harmonic. Since this is the only known 
analytic solution for the shape of a drag-reducing body at supersonic 
speeds, it is used as a standard to determine the accuracy of the numeri- 
cal methods used in the general method of this paper and also to determine 
the accuracy of applying the supersonic area rule to the nonslender con- 
figuration considered in this paper. 


INTRODUCTION 


The experimental verification by Whitcomb (ref. l) of the area rule 
has opened up a new re alm of possibilities of reducing the wave drag of 
complete configurations at transonic and supersonic speeds. The theo- 
retical basis for the area rule lay undiscovered for some time in the 
Hayes drag formula (ref. 2), which contains the implication that by 
indenting a body, a saving in wave drag equal to the drag of the wing 
panels can be produced. Thus, in principle, the wave drag of a slender 
combination utilizing an infinite cylindrical body can be reduced to 
zero. Lomax and Heaslet (ref. 3) have used this Hayes drag formula to 
determine wing-body combinations of minimum wave drag. They start with 
a wing alone and then, using the invariance principle of Hayes, distrib- 
ute multipoles along the longitudinal axis that will cancel the drag 
producing pressures of the wing alone. The wing-body combination thus 
obtained has zero wave drag if all multipoles are used. If, to obtain 
practical shapes, only a finite number of multipoles are used, the wave 
drag is a minimum for this number. 

In the foregoing method the problems of actual body shape and of 
pressure distributions are secondary considerations in the mathematical 
solution of the minimum-wave -drag problem. A method for minimizing wave 
drag of wing-body combinations dealing directly with the body shape and 
the pressure distributions has been developed based on the interference 
theory of references k, 5> and 6. The method, presented herein, is 
applicable to combinations with bodies which are approximately cylindri- 
cal in the region of the wing. The principal restriction on the wing 
other than small thickness and sharp leading edges is that it possess 
a horizontal plane of symmetry or that its upper and lower surfaces be 
independent . It is the purpose of this paper to present in the first 
part a general theory of drag reduction for wings on quasi-cylindrical 
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bodies. In the second part the method is applied to the reduction in 
wave drag of a combination employing a swept wing with sonic edges. 


SYMBOLS 



All lengths are measured in body radii, 
a body radius 

Aji aspect ratio of exposed wing panels Joined together 

c local wing chord or chord of rectangular wing 

c r chord at wing-body Juncture 

Cj) wave-drag coefficient of combination based on exposed wing 

area 


C BS 

DbCc) 

^(C)^ 

D C 

D S 

%(d) 

%(i) 

(1) 

Dtffi+d) 

(2) 

%(i+d) 

%(C) 

%(W) 

dz 

dx 


reduction in drag coefficient based on exposed .wing area as 
a result of body distortion 

drag of body in combination 

drag of body due to 2nth Fourier harmonic 

drag of combination 

saving in drag of combination by body distortion 

drag of wing in combination with body due to body distortion 

drag of wing in combination with body due to interference 

1 

drag of wing in combination due to interference and distor- 
tion pressures propagated along characteristic body cones 

drag of wing in combination due to interference and distor- 
tion pressures falling behind characteristic body cones 

drag of wing panels on distorted body 

drag of wing panels acting as part of wing alone 

local slope of wing surface in streamwise direction 




k 
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f 2 n( x ) 


^2Hg ( x ) 


g2n( x ) 



t 



li2n(x) 


ko^x), 
ko 2 (x), 
ko 3 (x) , 

l 


m 

M 


n 

2n 

P 

P 2 n( x ) 


P m 

00 

P 

^B(cL) 

^B(i) 

P B(W) 

p W(d) 


velocity amplitude function of 2nth harmonic neglecting tip 
effects 

velocity amplitude function for a pair of opposing sources 
inclined at the Mach angle 

amplitude of "body distortion due to 2nth harmonic 


successive approximations to g^x) 


nohhomogeneous part of integral equation for g 2n (x) 
components of ho(x) 

upstream limit of hody distortion 
slope of wing leading edge 
free -stream Mach number 
normal to quasi -cylindrical hody 
index of Fourier harmonics 
local static pressure 

half -amplitude of the 2nth Fourier component of the pres- 
sure acting on the hody due to the wing alone neglecting 
tip effects 

free -stream static pressure 

p - p m 

pressure coefficient, — — a 

hody pressure coefficient due to hody distortions 
hody pressure coefficient due to interference 
hody pressure coefficient due to wing alone 
wing pressure coefficient due to hody distortions 
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5- 


P W(i) 

^(w) 

q 

r,0 

r B 

R.P. 


Si(x) 


s Q ( x ) 


S P 

S W 

U,V,W 

v r 

V 

Vi 

V W 

w 2n( x >*) 


x,y,z 


x i 


X 2 


wing pressure coefficient due to interference 
wing pressure coefficient due to wing alone 
free -stream dynamic pressure, g-pV 2 

polar coordinates in y,z plane; y = r cos 0, z = r sin 0 

radial distance to point on "body of combination 

real part of a complex function 

semispan of wing-tody combination 

local semispan of inner portion of wing plan form 
(sketch (c)) 

local semispan of outer portion of wing plan form 
(sketch (c)) 

area of one exposed wing panel 
area of two exposed wing panels 
maximum thickness of wing section 

perturbation velocities along x, y, and z axes, respectively 

velocity normal to quasi -cylindrical "body, positive outward 

velocity along radius vector, positive outward 

free -stream velocity 

volume of tody indentation 

volume of exposed wing 

influence function 

coordinate axes with origin at vertex of wing alone, 
x measured downstream, y laterally to the pilot's right, 
and z vertically upward 

intersection with tody of most forward tody characteristic 
touching exposed wing 

intersection with tody of most rearward tody characteristic 
touching exposed wing. 
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x«,y»,z f 


Y(x) 

Z(x) 

z l 

zu 

P 

r2n( x # s ) 
6 

P 

‘Pb 


<P r 


CPi 


9, 


W 


oblique coordinates; x* = x - my, y* = y - mx, z* = z*7l - m I 2 

r w 0 (i)as 

J o 

r iw 0 (i)di 

J o 

local thickness of lower wing surface measured from chord 
plane, z = 0 

local thickness of upper wing surface measured from chord 
plane, z = 0 


•7m 2 - 1 
J W 2n (x,r)dr 


semiapex angle of wing section 

function defining the plan-form distribution of wing thick- 
ness (eq. (2)) 

constant, A 0 = 4, A n = 2 if n > 0 

free -stream density 

potential of distorted, body alone 


tip. 


potential of complete combination 
interference potential, cp c - (qfy + qig) 

potential of wing alone 
An asterisk (*) distinguishes a quantity associated with the wing 


I - GENERAL THEORY 

Statement of Problem and Determination of Boundary Conditions 

In references 4, 5> and 6 the boundary-value problem of determining 
the pressure distribution acting on a combination composed of a wing and 
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an infinite, circular, cylindrical "body was solved. In the method, the 
flow field of the wing alone is pictured as causing a virt ual distortion 
of the body surface, and the pressure field due to such distortion is 
computed term by term for each Fourier harmonic of the distortion. ' Fhp 
pressure field on the body is computed with the help of a table of char- 
acteristic functions or influence coefficients which give the pressure 
coefficient at any point of the body due to a unit bump at the origin. 

By summing the pressure contributions of each bump on the body by numeri- 
cal integration, the entire pressure field is built up. An essentially 
different and more complicated method using axial distributions of 
sources and multipoles was used to obtain the pressure field off the body. 
However, it has recently been possible to extend the method for computing 
the pressures for points on the body to points off the body by expanding 
the tables of influence coefficients. The theory of the new: me thod is 
set forth in Appendix A. This method which has been applied to combina- 
tions with cylindrical bodies ^having virtual distortions ■ is equally 
applicable to combinations having bodies with real distortions. Thus we 
can evaluate the pressure field and the pressure drag of a combination 
with a body of arbitrary small distortion using the above methods. The 
pressure drag of the combination can be minimized with respect to body 
shape by the calculus of/variations to determine the best shape for mini- 
mum drag. In this part of the report the theory of deter mining this 
shape will be considered. 

In this paper the wing alone is taken to be the wing panels ping the 
portion of the wing blanketed by the body when the leading- and trailing- 
edge intersections with the body are joined across the body in any arbi- 
trary manner. The wing -body combinations to which the analysis applies 
possess a vertical plane of symmetry and either a horizontal plane of 
symmetry or upper and lower wing surfaces that are independent (e.g., a 
triangular wing with supersonic leading edges) . An additional assumption 
is that the body shape at any point does not greatly deviate from a cir- 
cular cylinder ; that is, the body is quasi -cylindrical. Sketch (a) 
illustrates the winged portion of a 
wing-body combination together with the 
system of axes used in the analysis. 

All distances are based on a radius of 
unity so that the formulas in the anal- 
ysis can be written in simple form. 

The body cross section at any x posi- 
tion is taken as a circle on which are 
superimposed distortions that vary as 
cosines of even multiples of 0 to 
preserve a vertical plane of symmetry. 

Thus 

sL = (£)!>«=- ** « 

n=o 



Sketch (a) 
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The function ga n (x) is the shape amplitude function of the 2nth Fourier 
harmonic and is to he so determined that the -wave drag is a minimum. 
With regard to the wing, no limitations are made on its plan form. The 
coordinates of the wing section at any s panwise station y are taken to 
he 


zu = -z z = tmTiCx - pr + p,r) 


( 2 ) 


where T| is a function giving the thickness distribution and is zero off 
the wing. 

Consider now the "boundary -value problem of determining the combina- 
tion potential. This problem is solved in essentially the same manner as 
in reference 4. For a cylindrical body the combination potential is the 
sum of the wing -alone and interference potentials, but when the body is 
distorted, an extra potential q)g must be added. 


cp c = ^ + 9i + cPb 


(3) 


All potentials are to satisfy the wave equation in cylindrical coordinates. 


^9,1^, 1 a 2 (p d 2 q) _ Q 
dr 2 T dr ^ de 2 dx 2 

The first boundary condition is that the derivative of the combination 
potential normal, to both wing and body be zero with the body taken as 
the r = 1 cylinder and the wing as the z = 0 plane. Thus 


dcpp r = 1 

= 0 

dn r > 1 


0 < 0-< 2jt (body) 
e = 0,3t (wing). 


(4a) 


The second boundary condition is that there be no upstream waves. 


cpQ = 0 x < 0 


(4b) 


The second boundary condition is satisfied naturally by using Laplace 
transform methods, but the first one is not so easily met. Consider 
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sketch (t>) which shows a region of 
the quasi-cylindrical surface 
R(x,r,0) = constant in the neigh- 
borhood of point P. The direction 
cosines of the normal, to the surface 
at point P are 

Sr/Sx Sr /S r 1 Sr/Ss 

} , 

| grad R| I grad R| r |grad R| 

The velocity in the normal direction 
is thus 


z 



X 


Sr/Sx SR/Sr 5( Pc SR/Se 

| grad R | Sx jgrad R| Sr r2|gradR| Se 


For a quasi-cylindrical body the equation of the body surface is 


\ 00 x ' 

R(x,r,0) = r - Qy £cos 2n0 J g 2n U)<l5 = 1 


n=o 


so that 


Sr _ . Am'N Sr _ . Sr /W\ 

S^ " \ c v Sr " ’ Se W 

Since 

ox “ V ’ 3r ° \ G r/ Se 0 \?tJ 

the third term in the equation for v n is negligible compared to the 
first two terms so that 


Stp c Sq) c £r 


SR 


Sr 


Sx Sx “ " V Sx J r “ 1 


( 5 ) 
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Combining equations (l) and (5) yields the boundary condition for zero 
normal velocity at the body surface 


d<p c 

Sr 


dr dr dr \ c r/ Z_i 

n=o 


g 2 n (x)cos 2 n 0 ; 


r = 1 , 0 < 0 < 2 it 

(6) 


For no 


flow normal to the wing surface, we obtain 


dcp G /dz _ Szu 
V dx 


In terms of potential functions this boundary condition becomes 


1 ^C 
Sz r S0 


lf^w S( pi 


de + de + 




r>l; 0 = 0 ,* 


(7) 


Equation ( 7 ) is fulfilled by expanding cpj_ and cpL. in cosine series of 
even multiples of 0. Equation ( 6 ) is fulfilled by satisfying the 
following equations for each term of the Fourier series as in refer- 
ence 4. 


Sr 



n=o 


f Pri (x)cob 2n0 


^Pb 

Sr 


Y (N) X ^^ cos 2110 

n=o 


It is to be remembered that cp^ is the interference potential defined 
so that the combination possesses a cylindrical fuselage. It changes 
if the definition of the wing alone changes. All effects of distorting 
the circular cylinder are contained in cpg. 

Pressure Field Acting on Wing -Body Combination 


The first step in the present method for evaluating the drag of a 
wing-body combination is to determine the pressure distribution acting 
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on "the combination. This can "be done for the body very simply by rngyma 
of the W 2n (x) functions tabulated, in reference 5. The pressure field, 
off the body has been obtained, in references 4 and. 5 by distributing 
multipoles along the body axis. However, the determination of the pres- 
sure field, off the body can be reduced, to the same amount of work as 
determining that on the body by generalizing the W 2n (x) functions to 
W 2n (x,r) functions. This new procedure besides vastly simplifying the 
calculation of the pressures off the body also avoids the necessity of 
integrating across logarithmic singularities involved in the original 
method. The generalization of the W 2n (x) functions to W 2n (x,r) func- 
tions and their use in calculating the pressure field off the body are 
described in Appendix A. The W 2n (x,r) functions are plotted for 
n = 0 and 1 in figure 1. 

The pressure coefficients of the combination are the sums of the 
coefficients due to the body alone, wing alone, and interference. Accord- 
ing to Appendix A the pressure field due to the body alo ng (distorted) is 
given on the body and the wing as 


(Jgj ■ l X [«» w - * C 0“]~ ^ 


■Pff(d) 

(tm/cr) 



gg n (x- gr+p) 



x-p(r-i) 



-r 



( 9 ) 


The pressure field due to the wing alone can be determined as a pure wing 
problem in linear theory. Let 1365 the component of this field act- 
ing on the wing and the component acting on the body. Let *Ww) 

be composed of a component; due to infinite aspect ratio plus a component; 
due to the difference between infinite and the actual aspect ratio, that 
is, a component due to the tip. Performing a Fourier analysis of the 
■^B(W) field around idle body section and indicating the part due to the 
wing tip with an asterisk, we obtain 


P B(W) 

(iiu/c r ) 



n=o 


+ ?2n 



cos 2n0 


( 10 ) 
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The factor of 2 has "been introduced purely for later convenience, and the 
values of p 2n (x) and P2 n *(x) thus represent only half the true Fourier 
amplitudes. The quantity p 2n *(x) can "be thought of as due to the image 
system of line sources and sinks representing the tip ‘when the wing pres- 
sure field is constructed of line pressure sources in accordance with 
reference 7* 

The pressure due to the interference field is obtained in the same 
manner as in references !}• and 5, or analogously as that due to body dis- ' 
tortion. Let the normal velocity induced at the body surface by the 
wing-alone field be Fourier analyzed around the body in a manner similar 
to the wing-alone pressure field so that • 


Vr 

V 



n=o 



cos 2n0 


(11) 


Equation (ll) gives the virtual distortion of a cylindrical body neces- 
sary to conform to the wing -alone flow field. Since the cylindrical body 
is , in fact, not distorted by the wing flow field, it sets up an inter- 
ference pressure field equal and opposite to that associated with the 
virtual distortion. This opposite effect also explains the use of a 
minus sign in equation ( 11 ). Analogous to equations (8) and (9), we 
obtain 


*6(1) 

(tm/cr) 


p’ ^ ~j^an( x ) + ^2n*( x ) “ ^ J' jf2n(£) + 

n=o 0 


^2n*(^) J 




2nd 


( 12 ) 


j|V(l) : _2 y 1 J ^2n(x - fir + g) f 2 n*(x - pr + p) 

(%/c r ) P x/r + *s/F 



fsnU) + f an *( 6 ) W 2 n (j-r + l-|,r)dlj. (13) 
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Drag of Body and. Wing in Combination 


The resultant pressure acting on the "body of the combination is the 
sum of the pressure coefficients given hy equations (8) , (10) , and. (12) . 
Let the "body be distorted over the interval 0 < x < 2 . Then the body in 
combination has the drag 


A(c) ~I 0 l + + ^C 1 )] r E 


d0 dx (14) 


With r taken as unity, we obtain 


* t fo^r* i t 62 ”^ 08 h z { fen(x) + Bp - w 


P p 2n*( x ) + f sn(x) + tsn*(x) ~~f [g^d) + f snU) + 

f 2 n*(t) J W 2n ^Szl, 1^ dlj- cos 2n0 d0 


The total drag of the body in combination can be expressed as the sum of 
the drags of the various Fourier components since the components are 
orthogonal; that is, the pressures due to one component cause no drag by 
acting on the shapes of the other components. 


(dAO-q(c) ^ ( D /g) B (c) gn 


“I 


(tm/cy) Z_j (tm/c r ) 
n=o 


Then for n = 0 


(tm/ c r) 


= J + PPo( x ) + ft?o*( x ) + f 0 M + *o*(x) - 

3 f X [go(l) + fo(l) + f Q *d)] w o l)d6 j- dx (17) 
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and for n ^ 0 


(D/q), 


J S 2n (x) -{g^x) + Pp 2n (x) + PP 2n *(x) + f 2n (x) 


WM - | / [ffen(» + tsa.lt) + f 2n *(5)]w2n ( £ j i ^ dS } to 


The pressure coefficient acting on the wing in combination is given 
hy the sum of P w ( d ), Fft (yj), and Pw(i) • Consider the drag of the wing 
p ane l I>w(w) known from wing-alone theory, and evaluate %(&) and %(!) 
as follows: Consider first the pressure components 2g 2n (x- pr+p )/P Vr" 

and 2f 2n (x - pr + p) /p *fr which represent Ackeret pressures or two- 
dimensional linear-theory pressures attenuating inversely as the square 
root of r along the body characteristics. Including tip effects and 
considering both sides of both wing panels , we obtain for this component 
of the wing wave drag 


(p/q) w(i4d) = 8 r r y ^(x-pr+p) 

(W c r) P ^ L Vr 


3 P n=o 


.(x-pr+p) 


^ P / q ^W(i+d) 

(tm/c r ) 2 


F an*( x -Pr+P) 

*sF 


■ifJSJ 


ganCx-Pr+p) f. 

•7F + 


2n (x-pr+p) 


F sn*(x-pr+p) 

*«/r 


STi(x-pr+p,r) 


dr d(x-pr+p) 



D 
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It is convenient to integrate first out- 
ward along strips in the characteristic 
direction and then to sum "the strips 
between the leading and trailing char- 
acteristics tangent to the wing plan 
form as shown in sketch (c) . 


( D / q )w(i-KL) 

(tm/cr ) 2 




Sketch (c) 


fjaaCx-pr+p) + f2n*(x-pr+p) J J 


S ° dq(x-pr+p,r) 

Ct . d(x-pr+p) 


bx Jr 


( 21 ) 


or 


(D/qJ 


,(l) 


W(i+d) 8 


(tm/c r )‘ 


■ § 1 C [«“«♦'-« +f ^ w ]C cr ^5? ^ 


( 22 ) 


Consider now that part of the drag due to the pressure disturbances 
lying behind the Mach cones on which they originated and represented by 
the integral forms in equations (9) and (13) . This second component of 
the wing wave drag is 


(D/q) 


( 2 ) 

W(i+d) 




x-P(r-i) 



+ ^2 n(£) + 



(tm/c r ) 2 
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The volume of integration is shown in 
sketch (d) . In this sketch the posi- 
tive r direction has been reversed 
for clarity. Let us reverse the 
order of integration so that in the 
x,r,i space we integrate over a slab 
parallel to the x,r plane and then 
sum the slats with respect to 
Since W 2n (x,r) = 0 if x is nega- 
tive, there is obtained 


(Ww( i+ d) 

(tm/ c r) 2 


tip [-“> 


n=o 


+ 


f 2 n(0 + f 2n *(l)]^ // ^|- r *-|> r ) Cap ^ ar4 . (X “ PP+P) 

S P 

(24) 


or 


(h) 

(D/ q ) W ( 1+ d) 

(*fcm/ c r) 


00 _ Xo 

m 

n=o 


^S 2n ( 


(x) + f sn (x) + f 2 *00 to 


>] 
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The entire drag of the wing in the presence of the distorted body is 
thus 


( D /q) W (c) 

(tm/cr) 


( P / q) ¥(¥) 

(W c r ) 2 





dx 



dr](x,r) dr 
dx n/t 



+ f 2n (x) + f 2n *(x) I dx 







dr d(g - pr + p) 


(26) 


Body Shape for Minimum Drag and Drag Saving 


In this section we will derive an integral equation for the body 
shape for minimum wave drag. As noted in the Introduction, the alternate 
method of Lomax and Beaslet also yields a solution for body shape for 
mi ,ni mum wave drag. The two solutions are equal to the order of quasi - 
cylindrical theory. The solution by the alternate method as well as its 
relationship to the present method is described in Appendixes C and D. 

The drag of the combination can be expressed as the sum of the wing- 
alone drag plus the drag of a number of independent Fourier components. 

¥e can minimize the drag of any component or combination of components 
independently of the others; for instance, we can minimize the drag of the 
second component without altering the first component which controls the 
volume. The equation for the wave drag of the combination to be used in 
the minimization is 
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(D/q) p y~i r 

— - ^ / "ha. g^jiCx) -Ig^Cx) + Pp 2 n (x) + 3 p 2 I 1 *(x) + f2n( x ) + f 2n*( x ) ~ 

(%/cr) 2 P £i 4) 1 


J / X [652n(S) + f 2n (E) +f 211 *(S>] + ' 

IZJ" [fentx) + f 2n( x ) + f 2n*( x )]y °r + 

n=n ® ®i 


( D / q ^ W ( W ) 

(tm/c r ) 2 



+ f sn( x ) 




W 2n 


(|-r + l. 



dr d(g -pr + p) 


( 27 ) 


The range of the integrals in the terms was increased from x i < x < x 2 
to 0 < x < l with the help of simple properties of W 2 n (x,r) and T)(x,r). 
The constant >\ n is 4 for n = 0 and 2 for n > 0. In the minimization 
of the drag as given hy the above formula, the condition is considered in 
■which the shape of the indentation for the particular Fourier harmonic 
is at its value for minimum drag and then a variation 6g 2n (x) from this 
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shape is considered. The additional drag due to the variation, 8D C , is 
zero and is given as 


5(D/q) c g v, r 1 

(W=rf ' s A i 



+ 


J f X [%n(5> + W*> + ^n*(S)] l)«j- fa - 

00 2 ^ 

^ ^ g2n(x)dx ^Sgn ( ^ )^2n ^ dH + 

n=o ° ° 



^Tj(x,r) dr 
dx 

dx -Jr 



J jf W 2 n(^ - r + 1 - |,r)c r |l dr d(£ - |3r + p) ( 28 ) 

With the inversion of the following integral 

( 29 ) 
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equation ( 28 ) becomes 


^ pi f 

(Wc r f ** P L ^ J j^nM + PP 2 J X ) + 0P 2n *( x ) + W x ) + 

n=o ° 


- | jf g 2 n(5)W 3n |i^i|a{ - i J Jf 2I1 (|) + 


y- N 00 2 

w 2 n ^ J' 6g an (x) 


n=o 


C-^s-j// '- d— i-o 

Sp 


c r 


*3. 

81 


dr d(§ - pr + p) 


cbc 


( 30 ) 


Now 5(D/q) c must be zero to the first order and since 5 ^ n (x) is 
arbitrary and not zero, we must have for each n 
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p ^n-|2g 2n (x) + Pp 2n (x) + 3p 2n *(x) + f 2n (x) + f 2n *(x) J g gTt CljWpy, 


x-i 


ai- 


tf [w»+fEn*(»]^(^l)«}+| 


S ° dq (x,r) dr 


c r 


s i 


dx "Jr 


p J ' J w 2 n 0 “ r + 1 " p> c r g^-<*r ^ " Pc + P)1 =0 (31) 
Sp J 


This integral equation defining the optimum "body distortions for each 
Fourier component can "be written 


wherein 
h 2n (x) = 


g 2n (x) = hgnCx) f S 2n (£)W 2n 


x-| 


d£ 


(32) 


| [flPan W + PP 2n *(^) + f 2nW + * m *<*>] W ^n (| -*+ 


l-f,r)e r Mk^tg^ dr «*e-pr+p) - F* 0 Cr M^i 5^ 1 + ~ ' 

P ^ &5 J Si Sx ^ J 

2j3 ^ jj^an^) + f 2n*(^) J W 2n ^ p ^ ^ ^ (33) 
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The norihomogeneous term hg n (x) must "be evaluated for the particula r wing 
■under consideration, the values of p 2n (x), p 2n *(x), f P n (x ) , and f 2 „*(x) 
following from a kno wledge of the wing-alone flow field. Equation (32) 
is then solved "by numerical means to determine the "best distortion of the 
"body for the Fourier harmonic under consideration. The actual distortions 
are proportional to the wing thickness "because of the use of quasi - 
cylindrical theory and the application of the "boundary conditions on the 
r = 1 c ylinde r. The wing thickness must "be sufficiently thin, therefore, 
for the distortions to "be within the compass of quasi -cylindrical theory. 
For a wing of greater thickness an auxiliary condition can "be imposed on 
the allowable magnitude of the distortion and the minimum problem resolved 
subject to such a condition in a manner similar to that of Lom&x and 
Heaslet (ref. 3 ). 

A simple rule for the drag saved can be determined by evaluating the 
drag of the c ombinat ion at the minimum-drag condition. Under this condi- 
tion the drag is 


(D/q) c 

(tm/c r ) 


f X 1 S2n ^[" g2n ^ + p f S2u(^ W 2n + 

n=o 0 x 


(P/q) W ( W ) ^ 

(tm/cr) 


8 

P 


m 


f2n( x ) + -^! 


an 


*(x) 


n=o 


K' 


c r 


dq(x,r) dr _ 
8x ‘Jr 



r + 


1 . x , c r dr d(l - pr + p) .(3^) 

P / o| 

The last three terms represent the drag of the wing panels with a cylin- 
drical body, in which case the drag of the panels is the total drag of 
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the winged part of the combination. The first term thus represents the 
drag saved. _If the order of integration of this term is reversed, we 
obtain 


(D/q) s 

(tm/c r ) 2 


00 l 

^ ^n J S2n( x )|j32n( x ) “ ^ J 8en(£) w 2n^-jjpjI^cL!;l 

n=o 0 o J 


dx 



n=o 


*n 


f gsnWh sn 


(x)dx 


( 35 ) 


It is readily seen that the drag saving is simply the drag of the dis- 
torted body without any win g panels, and that this drag is the same 
whether the distorted body moves in one direction along its axis or in 
the other direction. If the combination is operated at a Mach number 
other than the design Mach number, the drag saving is less than that 
given by equation (35) • 


Body Shape for Minimum Drag by Alternate Method 


It is possible to build up the potential fields due to the addition 
of a quasi-eylindrical body to a wing, q>£ and cp B , by distributing sources 
and multipoles along the x axis as well as by the W 2n (x,r) method. 

Any set of solutions to the wave equation in cylindrical coordinates 
that is of the form of a function of x and r multiplied by cos 2n0 
can in principle be used for the purpose. In reference 4 a method is 
given of relating the body distortions to the multipole strengths for -the 
particular set of solutions used therein. The Bame method is used in 
Appendix B to show the relationship between body distortion and multipole 
strength for the multipoles of reference 3* In. Appendix C the equivalence 
between the quasi-cylindrical and multipole methods is used to develop an 
alternate method for determining the body shapes of a minimum-drag wing- 
body combination in a form amenable to analytical as well as numerical 
calculation. 


II - CALCUIATIVE EXAMPLE OF SWEPT-WING AND BODY COMBINATION 
WITH SONIC LEADING AND TRAILING EDGES 


An interest ing example of the general theory is the case of a swept- 
wing and body combination with sonic leading and trailing edges. By dis- 
torting the body properly it appears possible to reduce the drag rise 
that occurs as a swept -wing airplane approaches the Mach number for which 
its leading edge becomes 301110 and thereby to increase the Mach number 
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The supersonic area rule (refs. 8 and 9) gives an approximate idea of the 
shape of the distortion due to the axially symmetric harmonic. A compari- 
son is made in Appendix D of the "body shapes given "by the supersonic area 
rule and the present method. 


As an example, a minimum-drag computation is made for two harmonics 
of the swept wing and "body shown in sketch (e) . The calculation is for 
M = */2. If a combination at some other Mach number is desired, it can 
x he first transformed to M = •JW 

in accordance with Appendix E and 
then the answers should he trans- 
formed hack to the original Mach 
number. The results of the calcu- 
lation for swept wings with sonic 
edges are plotted in nondlmensional 
form applicable to any Mach number. 
A biconvex airfoil section has been 
used with a thickness distribution 
given by 



Sketch (e) 


- - (¥)£ ' ”) 


= tmTltx - (r-l) , r] 


(36) 


on the wing and 2 ^ = 0 off the plan form. The first major step in the 
calculation is to evaluate the various components of hyn(x) as given by 
equation (33) . The relative magnitudes of these components are in direct 
proportion to their importance in prescribing the proper shape for minimum 
drag so that considerable insight into the problem is obtained in this 
step. The second step is to solve the integral equation for the minimiz- 
ing shapes, and the third step is to evaluate the drag saving, volume of 
cutout, pressure distributions, etc. 


Velocity and Pressure Amplitude Functions 


The first four terms of equation 
normal velocity fields produced at the 



Sketch (f) 


(33) relate to the pressure and 
surface of the body by the wing- 
alone flow field. For the present 
example consider the wing alone to 
include the area blanketed by the 
body when the leading and trailing 
edges are extended as straight 
lines. The pressure and velocity 
fields of the wing alone can be 
built up by superimposing line 
pressure sources of the type given 
by R. T. Jones in reference 7. The 
u, v, and w fields of such a line 
source, as pictured in sketch (f) 
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are given in reference 10 for M = */2 as 


u 

y = 


me 


tVl-m 2 


cosh" 1 


J 


y* +z 1 



% — 7=== CORh " 1 r x> ' N )_ 

31 Wy 7a + z» 2 / 


cosh 


-1 


•VF 


+ 2 



f: cos " 1 

IT 


yy* +z 2 


>JyZ +z Z \ly1 Z + Z 1 Z 


wherein 


x* = x - my; y* = y - mx; 


s\/l-m 2 


For a sonic leading edge, let m approach unity and obtain the 
of equations (37), (38), and (39) as 

— - e_ ^x g -y g -z 2 
V " ’ * (x-y) 



Vx g -y g -z 2 

(x-y) 


cosh" 1 


*Jy z +z z 


2 : = 1 cos" 1 r _y_(y-. x ) + z2 ■ 

^ L (x- y)4y s +z 2 . 


(3T) 


r ] (38) 


(39) 


limits 

(40) 


(41) 


(te) 


From these equations the normal velocity and pressure distribution acting 
on the cylinder y 2 + z 2 = 1 can easily be determined. The velocity 
normal to the cylinder was computed for 15° intervals around the c ylinde r 
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and then Fourier analyzed using numerical methods for the pair of line 
sources 1-1 shown in sketch (g) as well as for a pair of tip sources 2 - 2 . 
The coefficients so calculated correspond to the Fourier series 


^ = f 0 jS (x) + f“ 2 S ( x ) c °s 20 + . . . 


C^3) 


v r * 

= f 0 g *(x) + f 2 s *(x)cos 20 + . . . 


m 


wherein the asterisk refers to the tip sources. The values of the 
Fourier coefficients are shown in figure 2. The constancy of the f 0 s (x) 

coefficient is noteworthy since 
the normal velocity distribu- 
tion varies with x. It is 
apparent that the tip has a 
minor influence. 

To obtain the velocity 
amplitude functions for the 
wing as a whole as defined by 
equation, (ll) , it is necessary 
to superimpose the effects of 
the sources and sinks that form 
the wing in accordance with 
sketch (h)'. The effect of the 
infinite number of small sinks 
between the leading and trail- 
ing edges is summed by inte- 
gration. The resultant values 
of f D (x), f 2 (x), f D *(x), p D (x), 
p 2 (x), and p Q *(x) are given in 
figure 3 . 



Mach cone intersections 

with body 


7 


-+++++- 



Infinite wing source-sink system Wing tip source-sink system 
Sketch (h) 


The main terms of equa- 
tion ( 33 ) are readily deter- 
mined. The final term involves 
numerical integration using the 
values of the W 2 n (x,r) func- 
tions in figure 1. The sixth 
term can be ascertained in 
closed form. Since from equa- 
tion ( 36 ) 
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we obtain 




2(x-l) 




The term involving a double integral can "be reduced to a single integral 
with the aid of equation (46) 

/ / W 2n (5 " r + 1 - x > r ) c rj (I - r + l,r)dr d(| - r + 1) 

J Dp 51 

= 2 JJ' W 2n (g - r + 1 - x,r) Tl - - ^J^l dr d(i - r + 1 ) 

Sp L -* • 

= 2 J jl . jj (^- r ) jd(g - r + 1) J W 2n (g - r + 1 - x,r)dr 


= 2 J (l ~ J w 2 n(l + 1 - x,r)dr 


With the following definition 


r 2 n( x > s) ./“ W Pr i (x,r)dr 
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we achieve the desired result 

rr we - r + 1 - v ) c Mk£iidai 4(8 - r + 1) 

^ ^Sp 

= 2 f (l - f-} r anU -x + ! r s)dl (50) 


Solutions of Integral Equations for First Two Harmonics 


The solutions of the integral equations represented hy equation (32) 
have "been obtained both by iteration and by simultaneous equations. If 

(o) 

a first approximation to g PT] (x) is known, say g 2n (x), then a second 
approximation can be determined from the integral equation as 


ffnvi ( ^ 


(51) 


The process is repeated until it converges (or diverges) . The conver- 
gence is greatly aided by a close first approximation. If a means of 
solving many simultaneous algebraic equations is at hand, the integral 
equation can be expressed as a number of simultaneous equations by writ- 
ing the integral approximately as a summation. The more equations, the 
more accurate the answer. The rate of convergence is decreased by 
increasing the interval 0 < x < 2 or by increasing the order of the 
harmonic. 

For the n = o harmonic the solution was started by using ho (x) as 
a first approximation. The component terms of ho(x) are shown in fig- 
ure 4. It is clear that the terms f 0 (x) and p Q (x) are equal and oppo- 
site and the tip effects are small, so that the effect of the body dis- 
tortions sending waves directly along the Mach cones is the dominant 
effect. This fact leads subsequently to a simple approximate theory. 

In the first calculation, the body was distorted only over the interval 
of the root chord. The subsequent approximations to go(x) were obtained 
by iteration in accordance with the equation 

f .\ _ n c+i (j) 

go (d+l) (x) =ho(x) +| J goU)W 0 |x-i|d| 


( 52 ) 
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The successive approximations converge so rapidly as shown "by figure 5 
that the initial approximation could, have "been used, as the basis of a 
good approximate solution. In the second calculation, the possibility 
of distorting the body in front of the wing was considered with the 
result that the calculated .distortion here is very close to zero, as 
shown by figure 6. For a general wing this corresponds to no distortion 
in front of the Mach rhombus enclosing the wing, a condition that follows 
readily from a consideration of momentum transfer through a large volume 
enclosing the combination. The effect of distorting the body downstream 
of the wing was investigated with the result that some distortion in this 
region is effective in further reducing the drag. However, the additional 
gain calculated for the present example is small. Thus, the significant 
interval for the purpose of the present example is the root chord, and no 
large change in the shape for minimum drag occurs over the interval of 
the root chord by considering a distortion interval greater than the root 
chord. 

A solution for the integral equation for n = 1 was attempted in the 
same way as that for n = 0, but it was soon apparent that the rate of 
convergence was relatively slow. For this reason a better initial approxi- 
mation was sought by solving a system of 11 simultaneous algebraic equa- 
tions using the trapezoidal rule to evaluate the integral. The approxi- 
mation so obtained is shown as g|^(x) in figure 7* Several iterations 
then produced convergence. 


Body Shape, Drag Saving, and Pressure Distributions 


The solutions to the integral equations, go(x) and g 2 (x), give the 
shape of the body for two harmonics from the integral of equation (l) . 


r-n r x 

— = 1 + (tm/c) J go ( I ) d| + (tm/cjcos 26 J g 2 (£)d£ 


( 53 ) 


The shapes of the indentations for the n = 0 and n = 1 harmonics at the 
wing -body junctures are shown in figure 8 for unit wing thickness. It 
is to be noted that the n = 0 indentation is similar to the wing profile. 
However, because of the finite chord length, the indentation does not 
come up to the full body radius at the wing trailing edge, a fact in 
accordance with the supersonic area rule. This point will be further 
elucidated when we consider effects of plan form on the n = 0 indenta- 
tion. The n = 1 indentation is much larger than that for n = 0, an 
effect that tends to discourage the use of the higher harmonics. For 
the 5 -percent -thick wing postulated, the indentation due to the second 
harmonic is so large that the assumption of a quasi -cylindrical body is 
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being stretched. The body cross section given by equation (53) incor- 
porating both harmonics is shown in figure 9* The c h a n ge in shape behind 
the tr ail i ng edge does not contribute much to the drag reduction as will 
be discussed later, and we therefore have the alternatives of incorporat- 
ing the n = 0 afterbody, n = 1 afterbody, the afterbody using two har- 
monics, or a cylindrical afterbody as in figure 9* The percentage of the 
wing volume removed up to the trailing edge is 69 percent and up to the 
end of the Mach rhombus (x = 8) about 101 percent. Body volume is 
removed only by the n = 0 harmonic . 

The slope of the indentation at the wing leading edge can be deter- 
mined from the integral equation with the aid of the fact that gg n (x) 
must be nearly zero in front of the Mach rhombus (see Appendix C) . The 
discontinuities in gg n (x) must coincide with the known discontinuities 
of hsnCx) at x = 1 because the integral of equation (32) is continuous. 
Hence, 


g2 n (l + ) = h 2n (l + ) - h2n(l”) 


(5*0 


The only term contributing to the discontinuity is the next to the last 
term of equation (33). With the help of equation (Vf) , we obtain 



Panels 
I 767 


/3C„ 

(tm/c) 2 

Body 


S2n( 1+ ) “ ^ (>/b 


(55) 


(I ) Exposed panels of wing alone 



Panels Body Harmonics 


A+B 

A+B+C 


5.48 

0 

n=0 

5-67 

0 

n = 0,l 


(2) Wing-cylindrical-body combination 



Panels Body 


Panels Harmonics and 
and intervals of 
body distortion 


A+B+C+0 

A+B+G+D+ 

E 


3.60 

.79 

4.39 


-.49 

2.61 

2.12 . 

n =0,l<^<S 

"■'•“ft 44 


13) Wing-distorted-body combination 
Sketch (i) 


Since each harmonic contributes 
a finite amount, the indenta- 
tion for h- 11 the harmonics is 
initially radial in the 0=0 
and 9 = it planes. 

The net drag saving due 
to body distortion has been 
computed by evaluating the 
drag of the distorted body 
alone in accordance with equa- 
tion (35) and also from the 
pressure distributions acting 
on wing and body in combina- 
tion. The distributions of the 
drag along the wing panels due 
to two harmonics of interfer- 
ence and distortion, shown in 
figure 10, have been integrated 
to obtain the net panel drags. 


The dr ag s of the panels and body in combination for various conditions are 
s ummar ized in sketch (i) which is to be studied in conjunction with 
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figure 10. Part (l) of the sketch shows the drag with the panels as part 
of the wing alone . Part (2) of the sketch illustrates how- the addition 
of the cylindrical body has the effect of reducing the panel, drag in the 
proportion 7*^7 to 5.67. (it has been assumed that the interference drag 
for harmonics higher than the n = 1 harmonic is negligible.) Part (3) 
of the sketch shows the drag saving due to two harmonics of distortion. 
Incorporating the n = 0 distortion into the body reduces the drag of the 
combination by 22 percent (of 5.67) and the addition of the n = 1 harmonic 
has the effect of reducing the drag of the combination smother 42 percent 
for a total of 64 percent. It is to be noted that distorting the body 
has the effect of reducing the drag of the wing in the presence of the 
body while increasing the drsig of the body in the presence of the wing. 

The wing has negative wave drag when two harmonics of distortion are 
incorporated into the body. The net gain in drag is the result of reduc- 
ing the wing-panel drag about twice as much as the drag of the body is 
increased. It is interesting that the section drag coefficients given by 
figure 10 for the harmonics of interference and distortion attenuates 
along the characteristics approximately inversely as the square root of 
the radius, so that drag savings can be realized for aspect ratios greater 
than that of the present example. 

On the body the drag loading can be defined from equations (17) 
and (18) as 


te ^ qs »4(c) 


(tm / c ) 2 


V 1 ^nfen( x ) / , . 

L s w i g2n(x) 

n=o 


+ P 2n (x) +P 2n *(x) +f 2n (x) +f 2n *(x) - 


f ffen(t) + *2n(£) 


+ WU)] W 2n (x- 6)« j- 


(56) 


This quantity, plotted in figure 11, for the n = 0 and n = 1 harmonics, 
shows the distribution of drag along the body. This drag results from 
the pressure field arising from the body distortion and by distorting the 
body in an existing pressure field. It is noteworthy that extending the 
interval of distortion aft of the wing trailing edge causes only minor 
reductions in wave drag for the first harmonic. An approximate calcula- 
tion shows this also to be true for the second harmonic. 

The pressure distributions acting on wing and body both in combina- 
tion can be readily determined, the distribution in the wing -body juncture 
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"being given by 


P 

(tm/ c ) 






[feiW - J •g 2n (5)W 2n (x-5)'i|J + [p 2 n(x) + P 2n *(x)]^ (57) 


The first component represents the effect of interference "between wing 
panels and cylindrical body, the second component the effect of "body 
distortion, and the third component the effect of the wing-alone field. 

The pressure distribution at the wing-body juncture for a cylindrical 
body is shown in figure 12(a) . The wing-alone component of this pressure 
distribution has an infinity at the leading edge of the wing-body juncture. 
Because the body is a perfect reflection plane as far as this point is 
concerned, the pressure coefficient is the same as that at the apex of the 
wing alone, namely, 8/31. The finitd contributions of a finite number of 
harmonics of interference cannot overcome this infinity. However, we can 
easily establish the curve for the combination near the leading edge by 
fairing it into the known value. It is to be noted that the resultant 
pressure distribution has a favorable pressure gradient for the boundary 
layer. In figure 12(b) the pressure distribution for a distorted body 
has been determined by adding the pressures due to two harmonics of dis- 
tortion to those for the cylindrical body. The resulting distribution 
has an unfavorable gradient. The effects at the wing-body juncture are 
typical of the wing as a whole since the pressure distributions on the 
span are essentially those of the root chord with an inverse variation 
as the square root of the radius (fig. 13 ) • They are also typical of the 
body since they are simply multiplied by a cos 2n0 factor. For a 
5-percent-thick wing the complete pressure field at the wing -body juncture 
shows a pressure-coefficient rise from about -0.20 at the leading edge to 
about 0.25 at the trailing edge. This rise of Ap/q = 0.4-5 (Ap/p = O. 63 ) 
is less than the sudden pressure rise, Ap/p, of 0.9 necessary to separate 
a turbulent boundary layer as given by Bogdonoff and Kepler (ref . ll) . 

It is thus probable that turbulent separation will not occur. 


Approximate Solutions for Axially Symmetric Harmonic; 
Effect of Plan Form 


On the basis of the relative magnitude of various terms of ho(x) in 
the integral equation for g^x), equation (32), it is possible to obtain 
simple approximate solutions for the indentation shape and drag saving 
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for the n = 0 harmonic. The relative magnitude of the various terms 
depends primarily on the wing chord, or more specifically on c/|3r, so 
that a short-chord theory and a medium-chord theory can he distinguished. 
In the analysis we will confine the indentation to the root chord and 
let PAe *> e unity orgreater so that no tip effects occur. The Mach 
number need not be J2 although the analysis is carried out for this 
Mach number. Under these conditions the effects of interference and the 
wing-alone pressure field represented by f Q (x) and p Q (x) are equal and 
opposite for this particular case, and equations ( 32 ) and ( 33) can be 
rewritten 


Sq(x) = ho(x) + | J g 0 (5)w 0 |x-i|di 


ho(x) = h^x) + ho 2 (x) + ho 3 (x) (59) 

with 




dq(x,r) dr 
dx Jr 


(6o) 


1° 2 « = JJ W 0 (5-r+l -x,r)c aT|(S ~* g +1,r) ard(l - r + 1 ) 


bo s M - | J* fo(E)Wo(- (62) 


The components of ho(x) involving the W 2n (x,r) functions are concerned 
with departure of the pressures due to body indentations from the values 
that would be calculated by linear two -dimensional theory, and as such 
these components assume importance for medium to long chord lengths only. 

Short-chord theory .- By neglecting hQ 2 (x) and hQ S (x) compared to 
h c , 1 (x) , we obtain an approximate simple theory that represents the solu- 
tion for vanishing chord length. Under this assumption the pressures 
acting on both body and wing are simply the linear-theory two-dimensional 
pressures based on the body distortion attenuated inversely as the square 
root of the radius. This simple physical model contains the primary 
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features of the flow due to the n = 0 harmonic, and the medium-chord 
approximation can he considered a slight modification of the simple model 
to account for three -dimensional effects on the pressure field. On the 
basis of the short-chord approximation, equations ( 58 ) to ( 62 ) yield 

--ij <*) 


Utilizing the result of equation (47), we obtain 


tO r=i - @0* w - - (;) &)[ l - *"■- » w 

and 


Thus the shape of the indentation is given as 

(1 - js) - i - 1! (1 - Wb - 1 ), (x,i) ( 66 ) 

The shape of the indentation is similar to the wing section, but the 
volume of the indentation is usually less than the volume of the wing. 

In fact the ratio of indentation volume to wing volume is 


Vj[ *fs - 1 

= s - 1 


(67) 


vo lume ratio decreases as the span increases although the actual vol- 
ume removed increases. 
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The drag saving calculated using equation (35) is 


(D/q) s = 




Eq 2 (x) dx 


( 68 ) 


Expressed as a drag coefficient "based on the area of the wing panels, 
the drag saving is 


Cp s _ (PMf) s ^ _ 32 (Vi-l) 
(W c ) 2 (W c ) 2 3* (-'/s + l) 


The coefficient of drag saving thus increases as the span increases "but 
so does the drag of the wing alone. The effect of chord length will sub- 
sequently "be discussed. 

Medium-chord theory .- It was possible "by neglecting three-dimensional 
effects to obtain very simple results for the first harmonic. However, 
for longer chord length than those to which the simple theory applies, 
account must "be taken of the terms containing W 2n (x,r) functions. The 
value of ho-^x) remains as "before. In the evaluation of ho 2 (x) a simple 
assumption permits an approximate determination to "be made readily. This 
assumption, expressed in physical terms, states that the three-dimensional 
effects represented "by the W 0 (x,r) function themselves attenuate along 
the "body characteristic in accordance with a simple approximation similar 
to that for the two-dimensional effects. Mathematically the approximation 
can he expressed as 


Wo(x,r) 


" ¥ 0 (x,l) 


W o (0,r) 

W o (0,l) 


(70) 


The validity of the assumption can he measured hy how closely the quantity 
W o (x,r)W o (0,l)/Wo(0,r) is invariant to changes in r. A plot of this 
quantity for various values of r is given in figure lh to show that the 
approximation is a good one. The value of W 2n (0,r) obtained from Laplace 
transform theory is 


Wan(0,r) 


l6n 2 + 3 
8*>/r" 


l6n 2 - 1 
8r s/2 


(71) 
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so that 


„ , , w 0 (xji) f 3 i ^ 

W Q (x,r) - + — J 


(72) 


From equation (72) we obtain the result 


w*> = -f f r 3 1 NA.a.V a, (73) 

^o ^o '-'x 8w o (0,l) Vo- 1 / 2 r 3 / 2 / \ c / 


Integration yields 


llQoW = 


(«Vi-l)(3»y~s+l) 

it-vTs 


| |l - 2(X e ~ l) jY(e-x + l)-§Z (c-x + l)j- j 


1 < x < c + 1 (74) 

with the following simple numerical quantities given in figure 15 : 


Y(x) W 0 (l)d|; 


r = 1 


(75) 


Z(x) = ^ iW 0 (£)dS; 


r = 1 


(76) 


Likewise the value ho s (x) with the help of 


fo(l) 




c + 1 


(77) 


can easily he expressed in terms of the foregoing numerical quantities 


as 


W x ) — §{ [ x - T ( x - x) + 1 z ( x - 1 )} 


(78) 


With a good approximation to ho(x) including three-dimensional 
effects, we can now solve the integral equation for the body shape. In 
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evaluating the integral we will approximate go(x) by ho(x) and apply 
the -trapezoidal rule to the interval 1 < x < c +1 broken into n parts 
of length c/n. Then 

n 

go(xi) = ho(xi) + ^ ^(^)W 0 |xi-|j| (79) 

3=6 

■where 


I 


3 



^ = \> 3 = 0 ,n 

Pj = lj 0 < j < n 


> 




( 80 ) 


The example already solved exactly for the interval 1 < x < c + 1 is used 
in figure 1 6 to measure the accuracy of the approximate methods . For a 
chord-radius ratio of 3 the short-chord theory is already inaccurate hut 
the medium-chord theory is in good agreement with the exact solution. 

Some improvement in accuracy can he obtained by iterating equation ( 79 ) , 
but for a more exact solution a better integration formula is needed - 
one recognizing the slope discontinuity in W Q |x;L - £i| Xi = li. The 
drag saving evaluated by the trapezoidal rule is 

°°S = (v) X (81) 

s 3=0 


Effect of chord-radius and semispan -radius ratios .- The effect of 
the geometric variables on the body indentation, drag saving, and volume 
of indentation have been calculated using the medium-chord theory to 
study the limitations of the short -chord theory. The drag saving is 
shown in figure 17 as a function of s/a and c/pa. The saving increases 
as the wing span increases and increases slightly as the chord increases . 
The short-chord theory gives a good rough measure of the drag saving up 
to values of c/pa of the order of 4. The volume of the indentation as 
a fraction of the wing volume is given in figure 18. It is seen that 
short -chord theory gives accurate volumes only for very short chords. 

Some actual indentation shapes are given in figure 19 . The main feature 
of the shapes is that for the longer chord lengths the indentations do 
not return to the original radius. This result would be predicted by the 
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supersonic area rule. Those that do tend to return do not have nearly 
as much volume as the wing panels. Thus, as the design Mach number is 
increased from unity, the leading edge remaining sonic, the body volume 
to be removed decreases from the wing volume at M = 1 to a fraction of 
the wing volume at higher supersonic speeds. If distortion is made over 
the entire length of the Mach rhombus, the transonic area rule gives an 
indentation volume equal to the volume of the wing panels at M = 1 and 
the present method gives 101 percent of that volume at M = >J2. 


Body Shape by Alternate Method 


The body shape for a minimum-drag wing and quasi -cylindrical body 
combination can also be obtained by the alternate method developed in 
Appendix C. For the example of this section, the alternate method gives 
the first harmonic body shape in an analytic form that is accurate to the 
or-der of quasi- cylindrical theory. This result is significant because it 
represents the only known analytical expression for the m-tn-tmnm drag body 
shape of nonslender wing-body combinations. As such it represents a good 
standard for determining the accuracy of the numerical methods used in the 
general method of this paper. It also is useful as a standard to determine 
how well the supersonic area rule predicts the minimum drag body shape for 
nonslender configurations. 

The details of determining the shape of the axially symmetric body by 
the alternate method are presented in Appendix D. The body slopes obtained 
are shown in figure 20 together with the results of the quasi-cylindrical 
theory. The two slope distributions agree to the order of the quasi- 
cylindrical theory assumptions. The body distortions corresponding to fig- 
ure 20 are shown in figure 21 for the winged part of the body. These two 
distortions are essentially the same but they differ considerably from the 
supersonic-area-rule distortions. Within the framework of theory, this 
difference is to be expected since the supersonic area rule assumes slender 
configurations and the present model is not slender. However, the practi- 
cal. consequences of this theoretical result must be determined by experi- 
ment because of the known Inaccuracy of linear theory for sonic leading 
edges . 


CONCLUDING REMARKS 


The Nielsen-Pitts theory of wing-body interference has been applied 
to the problem of minimizing the wave drag of wing-body combinations by 
utilizing radial body distortion. The method, applicable to comb ina tions 
with quasi-cylindrical bodies and horizontal planes of symmetry, deals 
directly with pressure distributions and body shapes. An integral equa- 
tion is derived for the body Bhape for minimum drag for each Fourier 
component of the radial body distortion. To demonstrate the use of the 
method, it is applied to the wave-drag reduction of a k-5° swept-wing and 
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body combination utilizing a wing chord of three body radii and a span 
of six body radii. The leading and trailing edges are taken as sonic, 
and the wing section is biconvex. The numerical results for this example 
support the following statements: 

1. The axially symmetric harmonic of body distortion causes a 
reduction in wave drag of about 22 percent of the drag of the wing panels 
on the undistorted body while the second harmonic causes a further reduc- 
tion of 42 percent. 

2. The over-all wave-drag reduction is -the result of reducing the 
wing-panel drag by about twice as much as the body drag is increased. 

In fact, with two harmonics the wave drag of the panels is approximately 
zero. 


3. The axially symmetric harmonic yields results qualitatively in 
accordance with the supersonic area rule in that the body indentation 
does not return to the full body radius at the trailing edge of the wing- 
body juncture . The volume of indentation is only about two thirds of the 
wing volume for an interval of distortion of the root chord. 

4. The magnitude of the maximum indentation given by the axially 
symmetric harmonic is about 50 percent larger than that given by the 
supersonic area rule. 

5. To the order of quasi-cylindrical theory the first harmonic body 
shape is the same as that given analytically by an alternate method. 

6 . The maximum distortion of the cross section at the wing-body 
juncture due to the second harmonic is about five times that due to the 
axially symmetric harmonic although no volume is removed by the second 
harmonic. This limits the wing thickness for which the second harmonic 
can be fully utilized to about 5 percent or less . 

7. The interval over which the body distortion was effective in 
producing reduction in wave drag is that of the root chord although some 
further small saving could be obtained by distorting the afterbody. 

8 . The pressure distributions produced on the wing by body distor- 
tion, although favorable for the wave drag, are unfavorable for the 
boundary layer . 

9 . A very simple short-chord theory can be obtained for the axially 
symmetric harmonic by neglecting three-dimensional body effects . An 
approximate treatment of these effects yields a simple medium-chord theory. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Feb. 7, 1955 
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APPENDIX A 

PRESSURE FIELD DDE TO DISTORTED BODY AT ANY SUPERSONIC SPEED 


Let us determine the pressure field due to a "body distorted in 
accordance with equation (l) 


ll] _ = Or) X g 2n( x ) C0S 2110 • (A1 ) 

n=o 

We will consider M = '/2 for the moment. The potential q>g is to 
satisfy equation (Al) on the "body and to fulfill the wave equation. By 
analogy with equations (12) and ( 16 ) of reference k, 


q>e = V 



Gyn( s ) Kan ( sr ) 
sK 2n f ( s ) 


cos 2n0 


(A2) 


wherein 


Gan( s ) =I'[g 2 n(x)] (A3) 

and Kan(s) is the modified Bessel function of the second kind. For the 
case r = 1 it was possible in reference 4 to express the inversion of 
equation (A2) in terms of characteristic functions independent of the 
shape of the humps of the "body. Thus in terms of the characteristic 
functions 

Wan(x) = L-riB&l + l] (A4) 

LK 2n »(s) J 

it was 


P = - 


shown that on the "body 
00 


= 2 g 2n (x) - J g 2n (i)W 2n (x- |)di cos 2n0 


n=o 


(A5) 
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It is possible to determine the pressure at any point off the "body 
in a manner similar to that for points on the body. This is accomplished 
hy generalizing the characteristic functions to functions of two variables 


ws,r) s + 1 ' 1 

L Kan (s) 


In Laplace transform notation the pressure coefficient is 


_ o ( V n tr.\ ^en( sr ) „ „ 
“ 2 VcF7 / G2n(s ' , cos 2n9 
\ r / IW f (s) 


00 

_ _p V -s(r-x) , ,f s(r-i) Kg n (sr) 1 x “I 

' w 'L + w 


With the aid of the following equalities from Laplace transform theory 


L e 


-s(r-x) 


Gan(s) J » gan(x - r + l) 


T -1 [ Kan(sr) e -s(r-x) _ n t , , 

L — — — — + ' _i __ — = WonCx - r + 1 ,r) 

JW(s) ^ 


,-l„ War) . e -=Cr-i) r x-r« 

““'“’[ww Hr J m J 0 fen«)W*-r + l-S, 

the final result for the pressure field on or off the body is 


r)di (A10) 


= l -f ' i,r)«]coB *» 


(All) 
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For any Mach number the pressure coefficient is 


-?&)£ 


n=o L 


g gn (x-3r+p) 


n x-pr+p , \ 

' p J g2n(^) w 2n(^-r+l-^>r^di 


cos 2n0 


(A12) 


In this equation the functions gpn(x) sire taken to he the actual ampli- 
tudes of the body distortion and do not depend on Mach number. If the 
combination is distorted as in the Prandtl -Glauert transformation (see 
Appendix D), then the g pn (x) functions will change with Mach number. 

To make computations using equation (A12) it is necessary to evaluate 
the integral numerically using tabulated or plotted values of W 2n (x,r) . 
Plots of W Q (x,r) and W 2 (x,r) sufficient for the purposes of this report 
are included in figure 1. 
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APPENDIX B 

RELATIONSHIP BETWEEN SHAPE OF QUASI -CYLINDRICAL BODY 
AND STRENGTHS OF AXIAL MDLTIPOLE DISTRIBUTIONS 


The method used in this report deals exclusively -with the flow field 
external to the wing-tody combination and requires no continuation of the 
solution inside the tody. However, it is intuitively clear that any 
quasi -cylindrical flow can te built up by superimposing distributions of 
sources and multipoles along the axis of the body. For the particular 
set of sources and multipoles used in reference 4 a set of characteristic 
functions, the M 2 n(x) functions, were defined and tabulated to allow the 
axial distributions of multipole strengths to be determined for a given 
body shape. (The problem of finding the body shape from the multipole 
distribution is a direct, although laborious, mathematical process.) The 
same characteristic functions are applicable to the particular set of 
multipoles used in reference 3- These multipoles written in several 
forms as follows 


cos 2n0 cosh (2n cosh 


-1 


-) 


^x 2 - r 2 




= cos 2n0L - KanC r s ) ] 


(Bl) 


are distributed along the axis in strength As^x) . The multipoles are 
started along the body axis at x = 0, as shown in sketch (j) so that 

r./Wx) 


x 



Sketch (j) 
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the tody distortion starts at x = 1 for M = */¥. The multipoles are 
not started at negative values of x since Apn(x) would not then have 
a Laplace transform when x < 0. The Laplace transform, of the potential 
due to the multipoles is 

00 

%(s,r,0) = - ^ ^ 2 n (s)K 2n (sr)cos 2n0 (B2) 

n=o 


and the potential in the physical space is 


^2, «x cosh ^2n cosh" 1 

q^(x,r,0) = - — } (cos 2n0) / Ag n (|) d | ( B 3 ) 

r,=n U o n/(x- 


S) - r 2 


Let q^. he equal to cpg produced hy a distortion of the "body given hy 
equation (Al) . Then, according to equation (A2) we have 


% (x,r,e) = = V^L- 1 £ cos 2n 8 (b4) 

n=o 


Comparison of equations (B2) and (b 4) establishes the relationship 
between body shape and multipole strengths in the Laplace space as 


^an( s ) = ” sKsn* ( s )-^2n( s ) 

2-sCri-baJcr) 


(B5) 
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Form the Faltung integral corresponding to equation (B5) with the help 
of the relation 


l -1 [SsLlfi 1 = L -i [Ws)] . x > x 

L 8 -i 

= 0; x < 1 


Thus 


g 2 n( x ) = " 


2itV r (t m /c r ) J o 


pX-l 

J ww 


cosh(2n cosh“ 1 x- i) 
7(x- I) 2 -1 


dl 


(B6) 


Equation (b 6) gives a method of determining the shape of a quasi- 
cylinder corresponding to a given distribution of multipoles along the 
axis. The reverse problem of determining the multipole distributions 
for a given body shape is not as simple. This problem will arise when- 
ever a minimum-drag wing-body combination is operated away from the 
design point. Equation (B5) is written in a form to facilitate the solu- 
tion of the reverse problem 


e -AanC 8 ) - “23iV 


CW 5 ) 

se s K2n*(s) 


(B7) 


The occurrence of Kgn * (b) in the denominator makes the determination of 
the inverse transform difficult if at all possible in terms of known 
functions since the zeros of Kan * ( s ) iu the complex plane are generally 
unknown. However, if we define the following function by its Laplace 
transform. 


Men( x ) 


L’ 1 







(B8) 



it 6 


MCA TN 3722 


then the solution to the reverse problem can he -written do-wn 


Apri (x - l) = -2iN 



x-i 

g2nU)Jfen(x- t)d| 


(B9) 


Tabulated values of M^(x) , M 2 (x) , Hi(x) , and Me(x) for 0 < x < It are 
included in reference 5 for numerical -work. Equation (B9) thus allows a 
simple calculation of the multipole strengths for a given shape. The 
functions M P n (x) have logarithmic singularities at the origin. Thus 
whenever g 2n (x) is finite, so is A 2n (x - 1) , and if g 2 n( x ) ^ ias a singu- 
larity, so will AenCx - l) . 
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APPENDIX C 

ALTERNATE METHOD FOR DETERMINING BOUT SHAPE 
FOR MINIMUM-DRAG COMBINATION 


In this appendix we will consider the solution for the "body shape 
for minimum drag "by the method of reference 3, and its relationship to 
the present method. In reference 3, the potential function for a general 
class of minimum-drag combinations is found without recourse to solving 
a boundary- value problem as already described in the Introduction. If 
the stream tubes corresponding to the potential are known, any stream tube 
can be replaced by a body to form a minimum-drag wing-body combination. 

To carry out the precise determination of the body shape requires the 
simultaneous solution of a pair of nonlinear differential equations. 
Practically, the shape must be calculated by some approximate method. 
However, the potential function for the minimum-drag combination is accu- 
rate to the order of linear theory. In the quasi -cylindrical theory it 
is assumed that the flow is known for a c ylindr ical body and wing combina- 
tion by the use of the interference methods of references 4, 5, a nd 6. 

Then the shape perturbation from the circular cylinder necessary to pro- 
duce minimum drag is calculated harmonic by harmonic. In this method the 
boundary conditions are applied on the r = 1 cylinder. However, for thin 
enough wing panels the minimum-drag body will approach a cylindrical body 
as closely as desired. (For both solutions the wing boundary conditions 
are applied on the z = 0 plane instead of -on the wing surface.) Under 
these conditions we may therefore equate the potentials given by the two 
methods to the order of quasi-cylindrical theory. 

In the present paper, the combination potential is 


9 C = % + ^ + 9b (Cl) 

and in reference 3 the combination potential is 

<P C = + % Cos) 

where cp^ is the potential due to the canceling multipoles. We can thus 
say that to the order of quasi-cylindrical theory 


(C3) 


% + 9b = % 
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Since from equations (B2) and (a 4) tlie following relationships are true. 


CO 

‘~k I cos 2nd A pn ( s ) Kpn ( sr) 


n=o 


00 

L(q>B) = V ^ cos 2nd 


n=o 


G 2n (s)K 2I1 (sr) 
s ^2n* ( s ) 


00 

L(q>i) = v (fe) X cos 2110 


n=o 


F 2 n(s)K^ n (sr) 
s K 2 n* (s) 


we obtain the following relationship among the -variables 


&2n( s ) - “F2n(s) - 


2rir(-b m /c T ) 


A pn (s)sKpn t ( s ) 


(c4) 


The following relationships permit us to form the Faltung integral 
corresponding to equation (c4) . 

L -X rwtfl] . L -1 [Ws) ] . cosh(2n GQsb~ 1 x) . X>1 

L S J 


= 0 


x < 1 (C5) 


The inverse of equation (C4) is thus 


S2n( x ) = -f2n( x ) - 


2jtV(tm/cr) 


fv ( i) cosh[2n cosh-^x-E)] d| 
v b <J(x - I) 2 - 1 


(c 6) 


The first term can be thought of as the -virtual distortion of the cylinder 
by the wing-alone field, and the second term is the body distortion due 
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to the canceling multipoles. Hie multipole strengths are determined by 
the method of oblique planes as used in the supersonic area rule and 
described in references 3 and 8 . 

The value of the shape amplitude functions, Esxyi x ) » can, in prin- 
ciple ■, be determined analytically from equation (C 6 ). The determination 
involves first calculating the normal velocity field at the r = 1 cylin- 
der due to the wing alone and performing a Fourier analysis to obtain 
the fpn(x) functions . Then the multipole strength distributions for the 
wing alone Agm(g) are obtained by the method of oblique cuts which 
includes a Fourier analysis. Finally, to obtain the shape amplitude 
functions, the integrals of equation (C 6 ) must be evaluated. 
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APPENDIX D 

AN ANALYTICAL DETERMINATION OF MINIMUM-DRAG BODY SHAPE 

BY ALTERNATE METHOD 


In this appendix the alternate method of Appendix C will be applied 
to the swept -wing example used in part II of the text to obtain an ana- 
lytical expression for g Q (x) . This is possible because f Q (x) can be 
expressed analytically for this example. 

For the purpose of this calculation equation (C 6) can be put into a 
more convenient form by choosing the inverse Laplace transform 

L” 1 [K 0 *(s)] = L _1 [-K 1 (s)] = j x > 1 

Jx 2 - 1 

= 0 ; x < 1 (Dl) 


Then for n = 0 the inverse of equation (C^) is 

Soto = -f 0 to /' X " 1 #o'(l) ■ - - ' • - 6 — dl (D2) 

2rtV(tni/cr) o «7(x- g) 2 - 1 



where 

2jt 

A°(i) = (d 3 ) 

As shown in sketch (k), 0 is the 
polar coordinate defining the Unp 
of intersection of the Mach cone 
with apex at § and a tangent 
plane; S^Cg.e) is the projection 
in the (y,z) plane of the wing area 
intercepted by the tangent plane to 
the Mach cone. An analytic expres- 
sion for the interference term, 
f 0 (x ) , has already been obtained. 

The first step is to compute 
^(1,0). The same wing alone is 
used as for the quasi-cylindrical 
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theory computation. Only the right wing panel need be considered since 
- S^ R (|,it - 0) due to the vertical symmetry of the wing. The 
equation for the surface of the right wing panel is 


Zu - ~z. 


■ ^ H*) 6 - h*) * 


y > 0 


The equation for the intercept of the tangent plane with the z = 0 plane 
is 

X = | + y cos 0 (D 5 ) 

Then the thickness distribution along this intercept is 

2% = 4^ e )] [1 . - °°° 8)1 (B6) 

Por the winged part of the body (0 < £ < c) 

V E ’ e) “ T ==-- a - / 2z “ as 

*v 1 + COS 0 S x 


■ [ j 


I - y(l - COB 0) I" _ I - y(l - COS B) 


da = s/(dx) 2 + (dy) s 


Sw E (i^0) 


I) ' (t) © + 


(1 - cos 0) - i (1 - COS 0) 2 j- 


2 
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the limit, y 2 , is a function of 


0 (see sketch (k)) and is given by 


y 2 = a 


0 < 0 < cos 



y* = 


1 - COS 0 


cos" 3 ^ < 0 < it 


The next step is to find AqCi) using equation (D3) . Use of the 

symmetry of the wing and the fact that (g,0) = - 0) gives 

Ij R 

Ao(5) V(E,e)ae - -*■/ '\'(i,e)ae = ^ f s^'d^Jae 


(D9) 


This integration must he performed in two parts due to the discontinuity 
in y 2 at 0 = cos -1 ^ From equations (d 8) and (D9) 

^1 - 1 + .J^ .^ 2sS I g2 j- (DIO) 


Differentiation of equation (DIO) and use of (D2) gives 



l) 

— d| 
(Dll) 


This e lli ptic integral is tabulated in reference 12. The interference term 
can he expressed analytically as 



MCA TS 3722 


53 


= . i fi - 2(x - 1) 1 (m 2 ) 

(We) * L c J 

When the integration for the present model (c = s = 3) is performed 

g Q (x) jL 2(x - i)~ | _ _8_ I TT 

(Wc) 7 I- 3 j " 3 7 - 


SLTJL e + jgjLjgE ll( a 2 ,k) 

2 x + 1 

k 2 = (5 - x)(x - 1) 

(7 - x)(x + 1) 

2 _ (X ~ 1) 

" (X + 1) 

here K, E, and H(cc, 2 ,k) are complete elliptic integrals of the first, 
econd, and third kinds, respectively. 
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APPENDIX E 

PRANDTL-GLA.UERT TRANSFORMATION FOR TOSG-BODT COMBINATION 


The calculations for the example configuration have all "been made 
for a Mach number of */2, hut the results can he generalized to any Mach 
number hy means of the Prandtl -Glauert rule . Let the configuration at 
Mach number M he transformed into an equivalent M = */2 configuration 
hy dividing its longitudinal dimensions hy p and leaving its lateral 
dimension unchanged. Let the potential at corresponding points he the 
same (see fig. 22). The two potential fields sire solutions to the dif- 
ferential equations given in the figure. Since the potential is the same 
in the y and z dime nsions, ckp/By and ckp/dz are unaltered so that the 
slopes of al l surfaces in the streamwise direction are unaltered. Under 
these conditions the quantities at M and at M = */2 are related as 
follows: 

Cr M = 

Hi = H/F 

b M = s >/a 

% ■ 

T % = <*%/? 

% = i VT 

At corresponding points x,y, z at M and x/p,y,z at M = */2 we have 


> (El) 
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The calculated results at M = V2 can be generalized to apply to 
all combinations of the family - in our case all combinations with 
sonic edges and fixed body radius and combination span - by plotting 


l"( r B / a) 


versus x/pa for body indentations, pP/(tm/c) versus x/pa 


p(tm/c) 

with r/a as parameter for wing pressure distributions, PCp/Ctm/c) 
versus s/a or Pita with c/pa as parameter for drag, and Vi/Vw 
s/a or pAg with c/Pa as parameter for volume of cut-out. 


versus 
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Figure 1.- Charts of characteristic functions. 
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(b) W 2 (x,r) 


Figure 1'.- Concluded 
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a = I 


Figure 4.- Components of ho(x) 
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Figure 5 -~ Successive approximations in solution to integral equation 

for n = 0. 
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Figure 6.- Dependence of go(x) on interval of distortion. 
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Figure 7*- Successive approximations in solution to integral equation 

for n = 1. 



Figure 8.- Contribution of axially symmetric and second harmonics to 
body indentation at "wing juncture . 
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A+B Panels on cylindrical body, n = 0 harmonic 


A+B+C Panels on cylindrical body, n = 0 and n = l harmonics 


A+B+C+D — Panels on distorted body, n = 0 harmonic 


A+B+C+D+E Panels on distorted body, n=0 and n = l harmonics 


Figure 10.- Spanwise distribution of drag on wing panels in combination 

with various bodies. 
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(a) Cylindrical "body. 

Figure 12.- Pressure distribution at wing -body juncture 
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(a) Axially symmetric harmonic. 

Figure 13 .- Spanwise attenuation of wing pressure coefficients due to body distortion. 
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Ob) Second harmonic . 
Figure 13.- Concluded. 
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Figure 1 6 .- Comparison of indentation by approximate and exact theories. 
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Figure 19. - Effect of chord and span on body indentation for n = 0 

harmonic. 
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(a) M >= 1 / 2 " (b) M 

Figure 22.- Relationship between wing-body combination and equivalent combination at ,M = >J2 





